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Abstract 

In the domain-wall formulation of chiral fermion, the finite separation be- 
tween domain- walls (L s ) induces an effective quark mass (m e ff) which compli- 
cates the chiral limit. In this work, we study the size of the effective mass as 
the function of L s and the domain-wall height mo by calculating the smallest 
eigenvalue of the hermitian domain-wall Dirac operator in the topologically- 
nontrivial background fields. We find that, just like in the free case, m e s 
decreases exponentially in L s with a rate depending on m,Q. However, quan- 
tum fluctuations amplify the wall effects significantly. Our numerical result 
is consistent with a previous study of the effective mass from the Gell-Mann- 
Oakes-Renner relation. 
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Chiral symmetry and its explicit and/or spontaneous breakings are important aspects of 
strong interaction phenomenology. Chiral dynamics dominates the low-energy hadron struc- 
ture and interactions. The chiral phase transition at finite temperature has been sought after 
experimentally for a long time. In addition, the weak interaction probes couple directly to 
the chiral currents, and the matrix elements of which sensitively depend on the chiral prop- 
erties of hadron systems. On the theoretical frontier, however, massless fermions defy the 
naive nonperturbative treatments. Indeed, for more than two decades, finding an appro- 
priate fermion formulation has been one of the most difficult challenges in lattice quantum 
chromodynamics (QCD). In the last few years, Kaplan and Shamir's domain- wall construc- 
tion and Narayanan and Neuberger's overlap fermion formalism || have emerged as 
promising approaches to simulating massless quarks. In this paper, we aim to study the ef- 
fectiveness of the domain-wall approach which has already been used in a number of realistic 
numerical investigations 

Following previous studies, we adopt Shamir's version of the domain wall fermion for- 
mulation 0, in which the five dimensional Wilson fermion is first introduced. The finite 
fifth dimension with L s lattice sites extends from s — to s — L a — 1. Dirichlet boundary 
condition on the quark fields is applied to the four- dimensional slices at s — — 1 and s = L s . 
The fifth component of the gauge potential is identically zero and the other four compo- 
nents are the same at every s slice. The above construction in the gauge- field-free case yields 
nearly-massless modes when the negative Wilson mass m (the wall height) is taken between 
and 10. These chiral fermions appear as the surface modes at s = and s = L s — 1 with 
opposite chirality. For < mo < 2, one flavor Dirac fermion can be constructed as 



where P± = (1 ± 7s)/2 are the chiral projection operators and n labels four- dimensional 
lattice sites. For the finite wall separation (L s ^ oo), the chiral mode on the s = wall 
couples with the one with the opposite chirality on the s = L s — 1 wall in an exponentially 
small way. Because of this coupling, a finite residual fermion mass is produced. The goal 
of this paper is to understand the size and dependence of this induced fermion mass on mo, 
L s when realistic background gauge fields are introduced. 

In the absence of the gauge potentials, the effective mass can be defined in terms of the 
pole of the free Green's function. For a large L s , it has a simple analytical form |2|], 



One can also obtain m e fj by either diagonalizing the hermitian domain-wall Dirac operator 
Hdw = 75-Ps-D or D^D = H\, w @>i]> where D is the domain- wall Dirac operator and P s is 
the reflection along the fifth dimension. On a lattice with the periodic boundary condition, 
the lowest four-momentum of a fermion is zero, and the lowest eigenvalue of '-/^PsD is just 
m e g. To study the effective mass in an external gauge field, we have constructed a code to 
diagonalize H DW numerically. To test the code, we have computed m eff on an 8 4 lattice with 
L s = 8, 12, and 16, and the result is shown in Fig. [I] together with the exact answer from Eq. 
(1). An ideal massless fermion is obtained at mo = 1 where m e g vanishes identically. For 
m ^ 1, the induced quark mass decays exponentially in L s with a rate depending strongly 
on rriQ. The exponential decay slows down significantly as mo approaches and 2. For any 
given accuracy e and any L s , there is a window surrounding m = 1 in which m c g < e. 




(1) 




(2) 



2 




In the presence of a realistic gauge potential, the effective quark mass result from the 
finite wall separation may depend on how it is defined. Different definitions shall yield 
results consistent up to a factor of order unity. One approach is to exploit the explicit quark 
mass dependence in chiral Ward identities such as the Gell-Mann-Oakes-Renner (GMOR) 
relation as done in Ref. 0. Here we explore the effective mass in an alternative way. In 
continuum field theory, the Atiyah-Singer theorem || states that the Dirac operator has a 
zero eigenvalue in the presence of an external background with topological charge \Q\ = 1. 
The explicit form of the solution was found by 't Hooft in 1976 0. On the lattice, however, 
the notion of topological charge is ill defined: any gauge configuration can be continuously 
deformed into a null gauge field. Moreover, the discretization of an instanton field can 
introduce finite lattice-spacing effects lifting any exact zero eigenvalue. Therefore, a test of 
the Atiyah-Singer theorem on lattice is usually complicated with various lattice artifacts. 

There exists, however, a definition of lattice topology and fermion zero mode which 
largely avoids this complication. In the overlap formalism, the Dirac operator is constructed 
from the overlap of two many- fermion ground states 0|. According to their recipe, one 
starts from a four-dimensional Wilson-Dirac operator with a negative Wilson mass m and 
calculates its eigenvalues. For m small and positive, the number of positive eigenvalues 
is equal to that of negative ones. When m increases, a level might cross from positive 
to negative or vice versa. When this happens, the gauge field is regarded to have a net 
topological charge \Q\ = 1. Then the overlap determinant is exactly zero by construction. 
This definition of lattice topology and zero mode do depend on, for instance, the Wilson 
parameters r and mo- However, the zero eigenvalue is exact, independent of the lattice 
spacing a and volume V. 
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The domain-wall formulation can be regarded as an approximation to the overlap formal- 
ism ||. Indeed, in the limit of L s — > oo, one recovers the overlap formalism apart from some 
unimportant discretization effect in the fifth dimension. For a fixed gauge configuration and 
Wilson mass m , if the overlap determinant is zero, Hdw nas an exact zero mode in the 
limit of L s — > oo. 

In short, in a background gauge field, if the hermitian Wilson-Dirac operator has a 
level crossing, the gauge field is considered to have a nontrivial topology. Edwards, Heller, 
and Narayanan have done extensive studies of the topological properties of lattice gauge 
configurations in this way [|HJ. In a topologically- nontrivial background thus defined, the 
domain-wall Dirac operator with a finite L s has small eigenvalues, nonvanishing only because 
of the finite wall separation. In the remainder of this paper, we are mainly interested in 
such domain-wall effects on the fermion zero-mode. We define the smallest eigenvalues of 
the hermitian domain-wall operator H^w as the wall-induced effective fermion mass. 




FIG. 2. Smallest eigenvalues of the hermitian Wilson-Dirac operator in a smooth 
(anti)instanton field described in the text. 



As a first nontrivial example, we have shown in Figs. § and |3| the results in a smooth 
instanton field configuration on an 8 4 lattice. A similar study has been reported in Refs. 
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The size of the instanton p is 10 and the cutoff parameter r max is 3. The flow 
of the small eigenvalues of the hermitian Wilson-Dirac operator is shown in Fig. 0. A 
level crossing from positive to negative is seen at m-o = 0.3. Four separate crossings in the 
opposite direction happen near mo = 2.2. Only three of the crossings are plotted. In the 
overlap formalism, the Dirac operator has an exact zero eigenvalue in the region between 
the crossings. 

The lowest eigenvalue of the hermitian domain-wall Dirac operator i^ow — 75-Ps-D is 
shown in logarithmic scale in Fig. |3|. The overall profile of the eigenvalue as a function of 
mo is similar to the free case in Fig. |l[ For a fixed L s , the smallest eigenvalue occurs at 
around mo = 1.3, shifted upward from mo = 1. This shift reflects the renormalization of 



4 



m cf f 0.0001 



le-06 







O 


+ 


+ 









□ 





o 


□ 


+ 


DWF Ls=8 










o 


DWFLs=12 










□ 


DWFLs=16 




□ 


□ 

















FIG. 3. Effective quark mass induced by domain walls in the same smooth instanton field as 
studied in Fig. [2|. L s is the number of lattice sites in the fifth direction. 



the Wilson mass in the presence of the external gluon field. As m deviates from mo = 1.3, 
the domain- wall effects grow stronger. For a fixed mo, the effective fermion mass decreases 
exponentially as L s increases from 8 to 12 and 12 to 16. 

Our result is quantitatively consistent with the chiral condensate (ipip) calculation in Ref. 
0. For instance, at L s = 10 and mo = 1.2, we have m e fj ~7x 1CT 6 . We expect then (ipip) 
grows like 1/my as the explicit quark mass parameter rrif reduces to ~ 10~ 5 . For a fixed 
rrif = 5 x 1CT 4 and L s = 10, there is a window in m ~ [0.9, 1.6] in Fig. |3| where the induced 
quark mass m e fr is smaller than rrif. Hence we expect (tpip) is approximately independent 
of m there. At m = 0.75, the effective quark mass is about 10~ 3 with L s = 12. We expect 
the chiral condensate to have little sensitivity on rrif when rrif > 10~ 3 . All of the above are 
in accordance with those reported in Ref. M. 

A smooth instanton on the lattice is far from a typical equilibrium gauge configuration 
entering in the Feynman path integral. A more realistic study of the induced quark masses 
requires gauge configurations with quantum fluctuations fully included. In the following, we 
work on a set of Monte Carlo configurations generated on a four-dimensional 8 4 lattice and 
with (3 = 6. The physical volume is somewhat small, but we suspect that the domain-wall 
effects have a weak dependence on it. 

We pick a lattice configuration in which the lowest eigenvalue of the hermitian Wilson- 
Dirac operator crosses from the positive to the negative at mo near 0.9. In Fig.|], we have 
shown the flow of the lowest few eigenvalues as a function of m . Between m = 2.35 and 
2.45, four levels cross from the negative to positive region. Although not shown explicitly, six 
level crossings occur at around m = 4. As pointed out in Ref. ||10j| , the hermitian Wilson- 
Dirac operator is symmetric with respect to mo = 4 and hence the flow diagram has the 
same symmetry. According to the overlap fermion formalism, the Neuberger-Dirac operator 
in the above gauge background has one zero eigenvalue when mo is in the interval [0.9, 2.35], 
three zero eigenvalues in [2.45,4], three again in [4, 5.55], and finally one in [5.65, 7.1]. 

In Fig.[5], we have shown the smallest eigenvalue of the hermitian domain-wall Dirac 
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FIG. 4. Eigenvalues of the hermitian Wilson-Dirac operator as a function of the negative Wilson 
mass mo in a topological gauge configuration generated on an 8 4 lattice and at /3 = 6.0. 

operator H DW . Different symbols correspond to three different L s = 8 (pluses), 12 (circles), 
and 16 (squares). Over a large region of mo, the effective quark mass decreases exponentially 
in L s , as is clear from the approximate equal spacings between pluses, circles, and squares. 
The fastest decay occurs at m around 2.0, compared with 1.0 in the free case and 1.3 in the 
smooth instanton field. The wall effects become strong again near m = 2.7 beyond which 
we find four small eigenvalues (not shown). To our surprise, this transition point to the 
doubler region is at higher m compared with the prediction from the spectral flow of the 
hermitian Wilson-Dirac operator. This may indicate some subtle differences between the 
eigenvalues of the transfer matrix in the domain-wall formalism and those of the hermitian 
Wilson-Dirac operator in the overlap formalism. 

The important point about Fig.[5] is that the magnitude of the effective mass is much 
enhanced relative to the case of the smooth instanton configuration. With a best choice of 
mo near 2.0, the domain-wall Dirac operator has the smallest eigenvalue (3 ~ 4) x 10 -4 at 
L s = 16. For the same L s and with mo of 1.3, the smooth instanton configuration yields an 
eigenvalue ~ 10~ 8 . This dramatic increase of the effective mass comes from the ultraviolet 
fluctuations. As we will discuss further below, the ultraviolet fluctuations at strong coupling 
can cause great trouble for the domain-wall formalism. In Fig. |6], we have shown the smallest 
eigenvalue of H DW in another topological gauge configuration. The result is qualitatively 
similar to that in Fig. |5]. In particular, for L s = 16 and mo = 2.0, the induced quark mass 
is around 10~ 4 . 

Since the effective quark mass depends on particular gauge fields, it is useful to get an 
average over an ensemble of gauge configurations. For this purpose, we have generated a 
set of 150 configurations on a 8 4 lattice at (3 = 6.0. By studying the spectral flow of the 
hermitian Wilson-Dirac operators, we found 12 topologically-nontrivial configurations. The 
second column in Table I shows mo at which the first level crossing occurs. All the crossing 
points are entirely concentrated in the interval between 0.8 and 1, a fact consistent with a 
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FIG. 5. Effective quark mass induced by the domain walls in the same topological configuration 
as studied in Fig. 0. 



similar study in a slightly larger lattice 8 3 x 16 |T0fl . We proceed to calculate the lowest 
eigenvalue of the hermitian domain-wall Dirac operator with L s = 8, 12, and 16 at mo = 1.8. 
The result is shown as the remaining columns in Table 1. We find the average effective 
quark mass at L s = 16 is 8(4) x 10~ 4 . The dependence of the average effective mass in L s 
is consistent with the exponential within the error. 

Our result can be compared with a previous study of the wall-induced quark mass using 
the Gell-Mann-Oakes-Renner relation on a 16 3 x 32 lattice at the same value of (5 0. In the 
domain-wall formulation, the quenched chiral condensate is related to the pion susceptibility 
X, by 



W> = K + m c G ff MOR )^ + b 



(3) 



where b is a constant vanishing in the L s — > oo limit. The effective quark mass m e R can be 
extracted from m/-dependence of the condensate and Xn- For L s = 16 and mo = 1.8, an 
extrapolation to m/ = limit yields m^ OR = 0.012. Note that the smallest at which 
the data is taken is 0.01. 

To understand the physical significance of m§j^ OR , we recall the spectral version of the 
GMOR relation in the overlap formulation HI 



rn 



fX-n 



mi 



' \Q\ 

mjV 



1 2(1-Af) ' 
V V A? + m 2 f 



(4) 



where ±Aj are a conjugating pair of eigenvalues of (j 5 D ov ) 2 (D m = overlap Dirac operator). 
The first term signifies the contribution from the topological zero modes; for sufficiently 
large volume V and/or large m/, this topological charge term is negligible. The second term 
measures the chiral symmetry breaking effects on the conjugating pairs of eigenvalues. In 
the study quoted above J7|, the topological charge term is insignificant even at the smallest 
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FIG. 6. Same as Fig. ||, but for a different topologically nontrivial configuration. 



iTif = 0.01, and the extracted m§j undoubtly measures the explicit chiral-symmetry 
breaking effects in conjugating pairs of eigenvalues as induced by the domain walls. 

One can push the quenched calculation in Ref. |7| to the limit rrif = 0. In this case, both 
the chiral condensate (ipip) and pion susceptibility are dominated by the zero modes. The 
quark mass obtained from the ratio of the two observables is just the smallest eigenvalue of 
the Dirac operator in the topologically-nontrivial configurations. Therefore in the quenched 
chiral limit, the effective quark mass determined from the GMOR relation coincides with 
what we have considered in this paper. To be sure, the effective mass defined from the 
effects on the fermion zero modes is not the same as the one defined from the effects on the 
conjugating pairs of eigenvalues. Nonetheless, both definitions shall be consistent within a 
factor of order unity. In this spirit, our average effective quark mass 8(4) x 10 -4 is indeed 
in accordance with 1.2 x 10 -3 from Ref. f7j. 

To be completely sure about the consistency of the two approaches, further numerical 
studies are needed. For instance, along the line of study in Ref. one can attempt to 
subtract the zero-mode contribution to the chiral condensate and pion susceptibility, and 
then both quantities can be measured in the rrif — > limit. On the other hand, the present 
study can be repeated at a larger physical volume; a 16 3 x 32 lattice will be more suitable 
for comparison. 

We finally return to the central question of the domain-wall fermion formalism: How 
large an L s is needed for a practical simulation? Clearly, the results for the free field or 
artificially smooth gauge fields are of no help here. The answer depends on the size of 
quantum fluctuations. For large values of (3 (and perhaps small physical volumes as well), 
such as the case we have presented, one can work with L s = 12 or 16 and keep the induced 
quark mass under control. However, in going to smaller /3 and larger physical volume, level 
crossing happens continuously in a region of mo above some critical value [K| . From Figs. |5| 
and |5], we expect that if a level crossing occurs slightly before the mo where the domain-wall 
Dirac operator is defined, the wall-induced quark mass will be huge. [When the level crossing 
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TTJUffl Tjo = 16) 
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14 


0.916 


fi 777e-03 


1 745e-03 


5 014e-04 


38 


0.874 


8 663e-03 


2 822e-03 


1 019e-03 


42 


0.876 


6 104e-03 


1 694e-03 

X . V;JiC WtJ 


5 601e-04 


45 


0.897 


5 849e-03 


1 382e-03 


3 584e-04 


•58 


0.945 


9 188e-03 


2 607e-03 


8 037e-04 


82 


0.862 


6.356e-03 


1.718e-03 


5.298e-04 


92 


0.901 


1.002e-02 


3.619e-03 


1.471e-03 


98 


0.853 


5.543e-03 


1.432e-03 


4.189e-04 


101 


0.842 


5.948e-03 


2.395e-03 


1.449e-03 


129 


0.845 


4.525e-03 


1.094e-03 


3.066e-04 


133 


0.877 


7.554e-03 


2.276e-03 


7.677e-04 


147 


1.003 


1.252e-02 


4.091e-03 


1.447e-03 


m eS 




0.0074(22) 


0.0022(9) 


0.0008(4) 



TABLE I. Average effective mass for the topologically nontrivial configurations 

happens slightly above tuq, the chiral symmetry breaking effects in the conjugation pairs 
of eigenvalues is expected to be strong. Again this leads to a large effective quark mass.] 
Then the average effective quark mass can be strongly influenced by this type of accidental 
configurations depending on the frequency they occur. The level crossings at large itlq 
reflect strong quantum fluctuations at the scale of the lattice spacing flIQ| . Therefore, it is 
not surprising that at (3 = 5.7, one needs to have very large L s (30 to 40) to keep m e ff small; 
of course, this is the price that one has to pay to keep the physical volume large. 

To summarize, we have studied the induced quark mass resulted from the finite domain 
wall separation by diagonalizing the hermitian domain-wall Dirac operator in topologically 
nontrivial configurations. We find the quantum fluctuation strongly enhances the domain- 
wall effects. However, the effective mass does show an exponential decay as a function of 
L s . Our result on an 8 4 lattice with (3 = 6 is consistent with the effective fermion masses 
from the GMOR relation, although a detailed analysis shows that the two definitions of the 
effective mass are not the same. Finally, we comment on the size of L s needed in a practical 
Monte Carlo simulation. 
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APPENDIX A: EIGENVALUES OF FREE DOMAIN- WALL FERMION 



The domain-wall induced fermion mass in the free case was first calculated by Shamir 
using Green's function approach. Vranas stated in his paper || that he obtained the 
same result by diagonalizing the domain-wall Dirac operator without showing the actual 
calculation. An explicit derivation of m e g in the rrif = case was later provided by Neuberger 



14]. Here we show a complete derivation with the inclusion of nif. 

The domain-wall Dirac operator in the free field limit in momentum space can be written 



as: 



D = ip + M+P+ + M~P_ , 
where P± = |(1 ± 75) and p = j^smp^. Mass matrices M ± are defined as 



(Al) 



(M-) 5 



$s-l.s' 



Kp)5 s 
b(p)6 s 



m f 5 s , Ns 5 s > :1 



[s, s 



(A2) 



(b 2 + l -b 
-b b 2 + l -b 



rrifb 




The second to (N s 



where b(p) = 1 — m + — cosp M ) and rrif is the explicit fermion mass. Our goal is to 
calculate the smallest eigenvalue of the bilinear hermitian domain-wall Dirac operator 

DD^ =p 2 + M+M~P + + M~M + P_. (A3) 

Since M + M~ and M~M + have the same eigenvalue spectrum, it is sufficient to consider 

\ 

.\r.\r= . (A4) 

\mfb ■■■ —bb 2 + m 2 J 

1) -th row of the secular equation (M + M~) ss '^/ S / = \ 2 ^/ s , or, 

(b 2 + l)* a - + # s+1 ) = A 2 ^ s , (A5) 

is solved by Hf a = exp[±as] (s = 1, • • • N s ) , provided A and a satisfy 

- 2b cosh(a) + (b 2 + 1 - A 2 ) = 0. (A6) 

The first and the last rows of the secular equation can be satisfied by a linear combination 
of exponential solutions, ^> s = exp[— a(s — 1)] + Aexp[—a(N s — s)] , where A is a constant 
to be determined: 

(b 2 + 1 - A 2 )(l + Ae- a ^-^) - b(e~ a + Ae~ a ^~ 2 ^) + m f b(e~ a ^ N ^ + A) = , 
(b 2 + m 2 - X 2 )(e- a ^-^ + A) - b(e- a ( N °-V + Ae~ a ) + m f b(l + Ae~ a( - N ^) = . (A7) 

Using Eq. (|A6|) to eliminate A from the above, we get 

2cosh(a)(l + Ae- a{Ns - iy ) - {e~ a + Ae~ a{Ns - 2) ) + m f {e-^ N "- 1) + A) = 0, 

(2 cosh(a) - (1 - m))/2b)(e- a ^- l) + A) (A8) 

_( e -*(N 8 -2) +Ae -a^ + + Ae -a( Ns -l)^ = q 



10 



Eliminating A from Eq. (|A8|) and rearranging terms, we have 

e 2a - e a /b - m){l - e a /b) + 2m f e- aN °{e 2a - 1) 

+ m}e- a ^~ 2 \e a -1/6) - e~ a2N °(l - e a /b) = 0. (A9) 

In the N s = oo, nif = limit, e a = 1/b. Assuming e a = 1/b + 5 and keeping terms linear in 
5, we find 

5 ~ -(m f + b N °) 2 (l-b 2 )/b. (A10) 
Finally, substituting e a = 1/b + 5 into Eq. ( |A6|) , we get the eigenvalue, 

A 2 = b 2 + 1 - b(e a + e~ a ) ~ -6(1 - b 2 )5 = (1 - fr 2 ) 2 ^ + 6^) 2 > 
or A = (1 - b 2 ){m f + b Ns ), as quoted in Ref. 0. 
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